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^ ■ Abstract 

! 

, The relation between an isotropic and an anisotropic model in loop quantum 

I cosmology is discussed in detail, comparing the strict symmetry reduction with a 

O^' perturbative implementation of symmetry. While the latter cannot be done in a 

. canonical manner, it allows to consider the dynamics including the role of small non- 

symmetric degrees of freedom for the symmetric evolution. This serves as a model for 
^ ' the general situation of perturbative degrees of freedom in a background independent 

^ . quantization such as loop quantum gravity, and for the more complicated addition of 

H ' perturbative inhomogeneities. While being crucial for cosmological phenomenology, 

■ - - ' it is shown that perturbative non-symmetric degrees of freedom do not allow definitive 

conclusions for the singularity issue and in such a situation could even lead to wrong 
claims. 



1 Introduction 

By far the most common and extended studies in cosmological models are those admitting 
only finitely many or even one gravitational degree of freedom, the scale factor or radius of 
the universe a{t) and maybe anisotropy parameters. This is justified in a homogeneous or 
isotropic model, symmetries present in our universe at large scales. This is also mainly im- 
plemented in models of quantum gravity, such as in Wheeler-DeWitt quantum cosmology 
|X 2} and loop quantum cosmology IH E] ■ That degree of symmetry does not only allow 
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insights into phenomenological aspects of the early universe but, in anisotropic models, 
also displays the most crucial features of classical singularities. 

More realistic models can then be obtained by breaking some of the symmetries imposed 
as before. After isotropy the first step is anisotropy, which can often be implemented 
exactly, and then inhomogeneity as the crucial one. In the first case (anisotropy) there 
is still a finite number of degrees of freedom and usually they do not pose much new 
technical difficulties. Conceptually, however, they are essential in particular for a general 
understanding of the issue of singularities (see, e.g., [H] for a recent discussion in the context 
of quantum cosmology). It is then of interest to discuss anisotropic models at different levels 
and to see what the implications for singularities are. Complete quantizations have been 
done in [71 IHl El with the result that their loop quantizations are free of singularities as in 
isotropic cases (TUl HH • 

Alternatively, one can try to implement anisotropics as perturbations to a quantized 
isotropic model and see how correction terms change the dynamical behavior. In this 
work we develop such a point of view both to shed light on the singularity issue and 
as a model for perturbative inhomogeneities. Following [71 where anisotropic models have 
been quantized in the loop framework, we introduce suitable anisotropy parameters in that 
model and implement conditions for them to be small. (As a first step we take only one 
anisotropy parameter of two possible ones which means that there is still a single rotational 
symmetry axis, a so-called locally rotationally symmetric (LRS) situation.) This will be 
discussed in detail at the kinematical and dynamical levels and results in an evolution 
equation describing an anisotropic model perturbatively around an isotropic one. The 
formulation can also be used for comparing background independent theories with theories 
formulated on a background, because the isotropic model can be seen as a background for 
the perturbative anisotropy. 

In section 2 we describe the classical models, to be quantized in section 3 both exactly 
and in a perturbative manner. This illustrates how symmetric models can be related to a 
less symmetric theory at the quantum level. In section 4 we discuss the evolution equation 
in terms of the perturbation and its new features compared to the isotropic model. In 
section 5 we look at the singularity issue from the perturbative point of view and show 
that conclusions drawn in such a situation must be used much more carefully than in a 
non-perturbative quantization. 

2 Classical models 

For simplicity, we consider the two most simple models, which are the fiat isotropic model 
and the Bianchi I model with one additional rotational symmetry (LRS). Since loop quan- 
tizations are based on Ashtekar variables [T^ IT^ given by a connection Al^ and a densitized 
triad E'f we first introduce these variables. In models, they take special forms correspond- 
ing to the fact that they are required to be preserved by a symmetry transformation up to 
a gauge transformation. 
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An isotropic connection and densitized triad for a flat model can always be written as 

Aa'dx"" = cdx' (1) 

in Cartesian coordinates, where the only two remaining components are conjugate to each 
other, 

{c,P} = Y^G (3) 

with the Barbero-Immirzi parameter 7 ^| ^] and the gravitational constant G. The 
variable p contains information about spatial geometry such as the volume 

^so = \P\'^' (4) 

while c is proportional to extrinsic curvature. Both variables are related dynamically as 
dictated by the Hamiltonian constraint 

3 

^iso = -^^r^c' V¥\ + ^^matter(p) = . (5) 



Similarly, in the LRS Bianchi 1 model we have 

A'^dx^Ti = Ar'^dx + Ar'^dy + cT^dz (6) 
d_ 

dx°- 

with symplectic structure 

{A, va} = 4777^ , {c, Pc} = Sn-fG 



Et—r' = PAr'd^+PAr% + Pcr'd, (7) 



There is a residual gauge transformation p^ \—>- —pA which can be fixed by requiring > 
(see 13 El for more details). The volume is now given by 



Kniso = \JPa\Pc\ (9) 

and the Hamiltonian constraint by 

i^aniso = -(87rG')-S~'(^'v^+ 2ACV^) + H^.tter{pA,Pc) = . (10) 

To identify the anisotropy as a perturbation it is convenient to introduce new variables 

(A,c) = (c + e, c-2e) (11) 

{PA,Pc) = {p + Pe, p-2pe) (12) 

together with the inverse transformation 

(c,p) = {l{2A + c),l{2pA+Pc)) (13) 

{e,pe) = {l(A-c),l{pA-pc)) (14) 

3 



as well as the symplectic structure 



{c, p} = y7G, {£, pJ^yTG. (15) 

The volume, in this manner, is given by p up to terms of at least second order in p^, which 
is the motivation for choosing our particular set of anisotropy variables: 



Kniso = V|P'-3P?P-2P3| ^ |^|3/2 (^1 - ^1 + 0(p,Vp')) 



(16) 



This shows that we have to assume -C p for perturbations, i.e. the approximation will 
break down close to classical singularities of the isotropic type where p = 0. Later, we will 
also have to assume £ -C 1. For the Hamiltonian constraint we obtain 

H^i^o = -^^\P\~'^'i^'p' - ^'P' - ^'Pe + '^CPePe + 0{pI/P^)) + H^.uer{P,Pe) (17) 

up to terms of third order in the perturbation. 

In the vacuum case, we have Hamiltonian equations of motion 

C = {C, -f^aniso} OC -C^p + - CEp^ (18) 

P = {P: -f^aniso} OC Cp^ - Cpl + pSPe (19) 

S = {£, i^aniso} OC \c^Pe - \cpe (20) 

Pe = {P£,-f^aniso} OC |cppe- (21) 

Thus, the isotropic variables receive corrections in their equations of motion only at second 
order. Ignoring those corrections, we recover the classical solution c = and p — const 
for flat Minkowski space. With this approximate solution for the anisotropic system, 
the equations for anisotropies simplify to e ~ and Pe ~ ~\p^e. The anisotropy in the 
connection will thus remain constant while that in the triad grows linearly in time until the 
approximations break down. Since the constraint is quadratic in connection components, 
a feature shared by the full constraint, we are not required at this point to assume £ <^ 1. 
However, the equations of motion show that p^ would be large if e is not small, and we are 
led to this condition if the perturbative evolution is to be valid for some time. 

We will later, at the quantum level, mainly consider an asymptotic approximation of 
very large p where all corrections of higher order in Pe/p are ignored. This also simplifies 
the classical equations because the Hamiltonian constraint is then simply proportional to 
— (c^ — e^)y/\p\. Notice that the quadratic nature of the constraint in connection variables 
makes this approximation reasonable as we do not need to expand in e. The solution e = c 
to the constraint is thus consistent even though the anisotropy e in the connection is not 
small compared to c (unlike p^ which must be small compared to p). As already mentioned, 
and used later, it is only necessary that e is small compared to one, which is usually the 
case also for c in semiclassical regimes. 
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Using the large-p constraint we have equations of motion 



c2 



P 

e 

Pe 



M 



2c^ 





+ 



2v^ 



With e = c = const we have p 



\t with constants 



of integration po and p^. Ehminating t, we have the internal time evolution 

Ve{p) =P° - \{P-Pq) 



(22) 



of with respect to p. Notice that this evolution, unlike the coordinate time evolution, 
does not depend on the value of e. We will use this fact later as a test of the perturbative 
quantum evolution. 



3 Quantization 

Loop quantizations are based on holonomies, and also in models one uses exponentials of 
connection components as basic expressions in addition to densitized triad components. In 
the isotropic case, this leads to a basic algebra given by p and exp(z/ic/2) for all G M, 
represented on a non-separable Hilbert space 

"^^iso = -^^^(l^Bohr, dyUn) (23) 

of square integrable functions on the Bohr compactification of the real line ^|, with 
orthonormal basis 

(c|/i) = exp(i/ic/2) , /i G M (24) 

and basic operators acting as 

exp(i/i'c/2)|yu) = l/i + yu') (25) 

M = ll^lM (26) 

using the Planck length ip = VSirGh. Similarly, for our anisotropic model we have the 
Hilbert space 

"^aniso — "^iso ® "^iso 

(27) 

with orthonormal basis 

(v4,c|/i, I/) = exp(i(/ic + z/A)/2) , z/ G M (28) 
holonomy operators as before and 

P^l/i, Jy) = i74/^l/^, J^) , u) = l-filiy\n, v) . (29) 
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These states are not completely gauge invariant in view of the residual gauge transformation 
fi 1-^ — yU, which however will not be relevant for our purposes. One can always ensure gauge 
invariance by working only with states J2fi. ,y'^^J.,Af^^ ^) where ip^^y is symmetric in /x. The 
triad operators then give us directly the volume operators for both models by inserting 
them into the classical expressions. 

For mathematical constructions one often makes use of dense subspaces Cyljso and 
Cylaniso two Hilbert spaces given by the so-called cylindrical functions which are 

finite linear combinations of the basis states or i/), respectively. With their algebraic 
duals of distributional states (linear functionals on the cylindrical subspaces) one has two 
Gel'fand triples Cylj^o C Hiso C Cyl*^ and Cyl^niso ^ ^aniso C Cyl*niso (see also (TE!). 

3.1 Symmetry reduction 

Since we are interested in the relation between a symmetric model and a less symmetric 
one we first use the opportunity to demonstrate how a reduced model can be obtained 
within quantum theory (see also [T71 UHl E] ) . We start with the anisotropic model and, 
following general constructions JTj, define an isotropic state as a distribution which is 
supported only on isotropic connections, i.e. e = 0. To fulfill the definition, we must find 
an antilinear map 

a:Cyli,,^CylLi,„,|/i)H-.(/i| (30) 

embedding the space of isotropic cylindrical functions into the distributional dual of the 
anisotropic one such that 

f^(l^))[|P>^)] = (/^|P>^)U=c for all |p,r). 

On the right hand side, we use the restriction of a state |p, r) to the subspace A = c upon 
which it can be interpreted as an isotropic state. Expanding cr{\fi)) =: XIk a '^K,A(p)(/t, A|, 
we have 

^P.r(p) = $^a,,A(/i)5.,p5A,. = ai\fi))[\p,T)] = I e-^^^/2e^(''+-)^/2d/iH(c) = 

K,A 

and thus 

(/i| = a(|/i)) = ^ <5p,«+a(k, A| = ^(k, - k| . (31) 

At, A At 

The summation is done over all real numbers, which is well-defined as a distribution in 
Cyl* ■ . 

amso 

In fact, in this case the definition simply amounts to multiplying an isotropic state in 
the connection representation with a 5-distribution supported at e = 0, 

K K 

but this would be more complicated in more general systems and the full theory where the 
definitions of |T7] still apply. 
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Indeed, this gives (/i|p, r) = 5;x,p+r, i-e. the evaluation of the distribution is non-zero 
only if the averaged label equals the isotropic one, 

and {ix\h^^h!f'^-^ = (/i|, i.e. exp(ip£:) acts trivially with its dual action on symmetric dis- 
tributions which are thus supported only on isotropic connections e = 0. 

Here, we already used the dual action of anisotropic operators on isotropic distributional 
states, which in this example mapped a symmetric state to a symmetric one. This is not 
the case for arbitrary operators, and one cannot simply define all operators for the reduced 
model by the dual action of anisotropic operators 18j. However, one can do this for special 
operators which suffice to define the basic operators of the model. For fiux operators, one 
can easily see that -|-pc = 3p is the only one which maps an isotropic state to an 
isotropic state: 

pa{\^J)) = 1{^\{2Pa + Pc) = li^lM = <y{f>\p))- (32) 

Moreover, it agrees with the isotropic fiux operator (f^H|) defined in the isotropic model. For 
holonomies, one can act with arbitrary products of Ha and he on isotropic states since we 
already know that hAh~^ acts as the identity and the remaining factor simply amounts to 
an isotropic holonomy operator. Also, all holonomy operators form a closed algebra with 
the isotropic fiux operator p = |(2p^ + Pc): 

[h^ihi^\p]\fi,u) = -l^il{p + r)\^^ + p,ly + r) = -l^i^p + r)h^P^h^;)\p,u) . (33) 

This is similar for the anisotropy operator which we can define by := ^ipA —pc)- Here, 
we have the commutator 

[hfhP,p;\ = -liiUlp - r)hfhP . (34) 

However, p^ does not map an isotropic distribution to another such distribution. 

If we are looking for a reduction of the operator algebra to implement isotropy, we 
need to find distinguished dual actions making use of the form of symmetric states. One 
clearly distinguished operator is p because it is the only fiux operator mapping an isotropic 
distributional state to another such state. For holonomies, hAh~^ is distinguished because 
it fixes any isotropic distributional state. Beyond that, no other operators are intrinsically 
distinguished, but we do not yet have a useful algebra of operators for an isotropic model: 
while p is available, our distinguished holonomy corresponds to the anisotropy e rather than 
the conjugate of p. Nevertheless, by acting with the distinguished p as well as all holonomy 
operators, we can generate the isotropic representation. More precisely, we define 

Stab(a) := {6 E Aniso|Oa(Cyli,J C a{CY\J} (35) 

and 

Fix(a) := {6 e AnisolOa(V') = for all ^ G Cyli,J . (36) 
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The isotropic algebra is then 



•ISO 



Stab((T)/Fix(a) 



(37) 



understood in the following manner: Fix{a) contains only holonomy operators besides 
the identity, and thus forms a normal subgroup of the Abelian group obtained from the 
multiplicative operation in the holonomy algebra. We can thus take the factor group 
of holonomy operators modulo Fix{a) which thanks to the distributive law can also be 
equipped with the usual additive structure to form an algebra. Due to the relation (jSSI), 
also the flux operators in Stab(cr) commute with Fix(cr) and thus descend to the factor 
group, defining Aiso- Using that Stab((j) is generated as an algebra by p and all h^^h^j'' 
and Fix((T) by h^h'^ , a general element of ^iso is 



with Laurent polynomials 7r„, allowing also inverse powers of hc^ for each n which are 
non-zero only for finitely many n (we can order p to the left using ()33|1 and factor out a 
power of hAh~^ into Fix{a) such that only the holonomy he appears explicitly). Through 
this construction we derive an algebra isomorphic to the isotropic one. 

As the presence of factor spaces indicates, there is no canonical isotropic subalgebra 
of the anisotropic algebra without additional input. This is analogous to the classical 
situation where the splitting of anisotropic variables A and c into an isotropic average c 
and an anisotropy is not unique. Classically, a given p, as we have it distinguished at the 
quantum level, determines a form of e as a linear combination of A and c by requiring 
{p,e} = 0. Indeed, also the quantum analog of e, given by hAh~^, is distinguished as it 
commutes with p. However, a unique form for c and Ps then follows only after an additional 
choice, which for us was the fact that the volume should receive corrections only to second 
order. We use this here to define Pe = lipA—pc) as above and then require that (jH^ vanish, 
analogously to {c,Ps} = 0. This distinguishes h\hc as the isotropic holonomy operator in 
addition to the flux p which indeed form a subalgebra of the anisotropic operator algebra 
mapping isotropic states to isotropic states in a way isomorphic to the isotropic model. 

In this manner, we obtain a unique subalgebra of the basic anisotropic operator algebra 
which is isomorphic to the basic isotropic algebra. This also presents an alternative way 
to define a symmetric model without referring to symmetric distributional states at all: 
we use the classical relations between symmetric and non-symmetric variables to define a 
distinguished algebra derived from the non-symmetric holonomy-fiux algebra. Using the 
fact that the non-symmetric representation is cyclic on |0,0), i.e. the subspace ^aniso|0,0) 
is dense in Tianiso, we generate the representation of the reduced model by acting only with 
the subalgebra on the cyclic state. This defines the reduced state space, which is equipped 
with an inner product by requiring holonomy operators to be unitary. Upon completion, 
this inner product space defines the Hilbert space of the isotropic model and its quantum 
representation in agreement with the loop quantization of the classically reduced model. 
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Characteristic properties of the representation are then clearly inherited from that of the 
less symmetric system, which in some cases is also visible for the reduction from the full 
theory [HIIH]. 

3.2 Perturbations 

Dynamical information is obtained by quantizing the Hamiltonian constraints and asking 
physical states to be annihilated by the resulting operators. Because the constraints contain 
connection components, and connections are represented via holonomies which shift the 
labels of states as in (j2SI), quantum constraint equations become difference equations jTHl 
inn 17] for states \ip) = XI/^'^mI/^) 1"^) ~ 'Yli^iv''^tJ-Al^^^) triad representation. For 

large scales, e.g. I'm. the isotropic case, the difference equation is well-approximated by 
the Wheeler-DeWitt equation [20] while for small there are crucial differences implying 
the absence of singularities at the quantum level fl]. In order to see the relation to 
the Wheeler-DeWitt equation, one has to assume that the wave function does not vary 
rapidly on small scales, i.e. changes in /i of the order one, which allows one to Taylor 
expand the difference operators into differential ones. This limit, however, does not exist 
for arbitrary states because there is no operator for c on Tiiso C2]- Indeed, the Wheeler- 
DeWitt quantization is based on a representation on the usual Schrodinger Hilbert space 
TYs = L'^i^i dc) which is inequivalent to the representation on the Bohr Hilbert space used 
for the loop quantization. 

In this paper, we analyze the case of an anisotropic model where, however, difference 
operators are not expanded for both triad components but only for that corresponding to 
the anisotropy Pe- We thus assume that the wave function does not vary rapidly when 
changes and that in an approximate sense an operator for e exists. Although we need to 
assume Pe to be small, this will not push us in the regime where discreteness is relevant 
provided that p is large. The key additional assumption is that e is not too big which now 
becomes important at the quantum level, and indeed both anisotropy parameters, p^ and 
e, have to be small. (The situation is different from pure isotropy where small p usually 
implies large c.) When only small e are allowed, we do not probe the whole configuration 
space and do not see its compactness which underlies the Bohr Hilbert space used in 
homogeneous loop quantizations. As in the large volume limit, we therefore introduce the 
Schrodinger Hilbert space Tis for anisotropics, which we do by defining the perturbative 
Hilbert space 

^pcrt := ^iso ® 7^8 (38) 

and realizing its dense subset Cy\^^ ® Cylg as a subspace of the dual Cyl^jjjgQ. In the 
Schrodinger Hilbert space we have to choose a suitable dense set Cylg, which for our 
purposes will be the set of all functions which are products of a polynomial and a Gaussian. 
Indeed, Schrodinger states in Cylg can be interpreted as distributions on the Bohr Hilbert 
space [21] using the antilinear map 

tt: Cylg ^ CylLhr,^(^)[|p)] := / e^^^^de . 
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Dual actions of basic operators on those states are given by 



{Pe<m\p)] = vr(V^)[pt|p)] = i74|pe'P^^(e)de=i74/e^'''-z^V;(e)^ (39) 
= 7r(-lz74d^/de)[|p)] 

and 

(e*"%(^/'))[|p)] = ^[e-'^'\p)] = j e'^P-^^'^de = j e'^'T^^^^de = vr(e'""V')[|p)] • (40) 

The first equation shows that the momentum operator is just the derivative operator, as 
also on the Bohr Hilbert space, while the second shows that {e'^'^^7r{ijj)){e) = e^'^'^ijj{e) where 
on the Schrodinger Hilbert space we can now take the derivative with respect to u. By 
going to the dual action on the image of vr, we thus obtain a simple multiplication operator 
for e well-defined with domain 7r(Cyls). 

In addition, if we choose an appropriate state \E' G Cylg we have an embedding id ® 
\Ef: Cyljgo — >• Cylpgrt of isotropic cylindrical states to perturbative cylindrical states where 
anisotropies are small in mean value, but not eliminated exactly. Moreover, since we have 
to choose a state \E', which in the above notation is identified with the map C — > Cylg, 1 i— * 
\E', there is no unique canonical embedding. (Coherent states are natural candidates, 
but unlikely to be preserved dynamically.) In addition to the strict implementation of 
symmetries a we then have the perturbative implementation 7r,j, := {'k^7c)o [id^"^) which 
are both maps from Cyl^go to Cyl^^jg^: 




We have the two Gel'fand triples of the exact models on the left and right, with the canon- 
ical embedding a of the reduced model through distributions implementing symmetries 
precisely. In between, we have the perturbative cylindrical space which is connected to 
the anisotropic side through the antilinear map vr from Cylg, combined with the antilinear 
dualization ^ of isotropic states. 

In the perturbative sector, non-symmetric degrees of freedom are present but can be 
arranged to be unexcited at least for initial states of the evolution. In this sense, the 
situation is similar to the field theory example studied in [22] where non-symmetric field 
modes are included in their vacuum or other coherent states. In this case, the field theory 
vacuum is a natural candidate for the analog of Perturbative states can then be thought 
of as states where only the symmetric modes have been excited but non-symmetric ones 
have been left in their vacua. This is technically similar to our construction of strictly 
symmetric states at the end of the preceding subsection where we used only symmetric 
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operators in order to generate all states of the symmetric model out of the non-symmetric 
cyclic state |0,0). Indeed, for background independent systems in loop quantum gravity 
the analog of a ground state is the cyclic state where no geometry at all is excited. Thus, 
although physically the construction in |22] is closer to the perturbative situation here, 
technically it is analogous to the situation of the strict model. This illustrates crucial 
differences in the context of symmetries between quantum field theory on a background 
and a background independent quantization of gravity. 

We have already discussed the relation between anisotropic operators and isotropic 
ones for the exact implementation of symmetry. For the perturbative implementation, the 
same relation exists for operators which do not depend on anisotropics, but we can now 
also act with other operators as we already saw for e and in pn]) and ()39p. In general, 
an anisotropic operator O acting on Haniso has a dual action on * ® 7r(Cylpgj.t), but does 
not necessarily fix this subspace of Cyl*jjigQ. However, it does so perturbatively when we 
expand it as a sum of operators in the perturbative sector. Choosing, for definiteness, a 
Gaussian for \l/ we have perturbative states of the form 

where c = |(2A + c) and e = — c) are understood as functions of A and c. By the 
chain rule, we then have flux operators 

and 

= -^i^iy = -^74 (I - 1) ^ = (h4^ - 2p.) ^ . 

In composite expressions such as the volume operator one can then expand functions of 
flux operators in when they are expressed in the flux representation, and analogously 
for holonomy operators in the connection representation. 

This is most easily done if eigenvalues are already known. For instance, if an operator 
O has eigenstates u), we take the eigenvalues O^^i, and insert 

/i = p + A-f-^ip^ps = p + P , 

y = ly- A-f-Hp^p, = \y-P (41) 

with dimensionless P := Ap^/jip, which follows from the classical relation (jl2p or the 
calculation above. This yields a function 0{u, P) which we expand in the perturbation P, 

0{u,P) = J20ili'^)P'. (42) 

k 

Note that P itself need not be small by our assumptions, which would mean p^ <^ £p. 
We have, however, P <^ u such that each olg^{u) must drop off at least as z/"'^. For 
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any fixed k, the values Ol^^i^u), interpreted as eigenvalues, define an isotropic operator 



Aso* ~ J2u^isoi'^)W){'^\ ^^ch that Ol^^lu) := Ol^^{i')\i') . Thus, we obtain the expansion 



where the right hand side acts on Tipert- This approximation of operators on a selected 
class of states implementing small anisotropy is discussed in more detail in App. ^ which 
can be applied to any Hilbert space |z/) (g) Hs for fixed z/. We understand here always a 
finite truncation to a given order in k, and do not discuss convergence of the expansion 
in whatever sense. We will later use states \E' of the above form, assuming a < Eq <^ 1 
in order to expand also in e. This results in similar expansions for operators containing 
holonomies, such as the Hamiltonian constraint. 

3.3 Hamiltonian constraint 

The Hamiltonian constraint operator for the anisotropic model results from the expression 
given in after using A = ci = C2 and c = C3 (as in j7j, we ignore an ambiguity parameter 
^21 which is of the order one and multiplies connection components in the constraint): 




(43) 



k 



H 




ilA) cos(iA) 6^} 
(44) 



sin (|c) y cos (|c) — cos (|c) V^sin (ic) 



and similarly Oa by replacing c by A, where V is the anisotropic volume. Making use of 
the change of variables (fT^ one gets 



+ 2 sin (ic - e) cos (^c - e) sin (|c + ^e) cos (^c + ^e) (45) 

Here Oc and Oa are more conveniently dealt with by using their action on the u) 
basis, namely. 




Now following we expand Oa and Oc up to second order in p^, using from the exact 
analysis [Zj that 
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Inserting (jUJ, we obtain the eigenvalues 

(O^W,PWP,P) = ¥ (|74)'^' {\lv + 1 + P\-\lv -1 + P\) J\\v - P\ 



For P -C z/, the absolute values can be evaluated as 



\lv + l + P\-\lv-l + P\ = ||z^ + l|(l + (fP + l)-ip)-||z>-l|(l + (|z>-l)-ip) 

= Ifz/ + 1| - Ifp - 1| + (sgn(|z/ + 1) - sgn(|p - 1)) P 



2sgn3/2(z/) - 3X3/2 (i^)i^, 



introducing 



Xsi.^^) ■■= -^(sgn(z/ + 5)- sgn(z/ - 5)) = | ^^^^ ^ ^ (49) 



sgn,(z7) := l(|zy + 5|-|zy-5|) = <( r^z? for|z7|<5 , (50) 



1 ioi v > 5 
for |z/| < c 
-1 for v < 5 



and 



||z7 + P| = ||z7|(l + 3P/2z7) 
while the other expressions are expanded as 



= yiz^ (1 - |P/z7 - fPVz/^ + 0(PVz>3)) 

and 



1^7 + l-P|V2_|l^_l_p|l/2 = |lp + l|l/2 /l_(lp+l)-lp 



-|ip-i|vyi-(iz7-i)-ip 

+ 1|1/2(1 _ + l)-lp _ + l)-2p2) 

~\\v-i\ni-w-^)-'p-w-i)-'p') 

+0(PVz/3) 

|l^,+ l|l/2_|l^_l|l/2 

-|(||z>+ iri/2gg^(iz>+ 1) - \\u-l\-^'\gn{\u-l))P 

+ l|-3/2 _ _ l|-l/2)p2 + 0(PVz>3) 

2V3A3V^-3v^PA3^^ 



-|v^p2A3|z7|-3/2 + 0(PV 
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introducing, more generally than before, 

Asf{u):=^^{f{u + 6)-f{u-6)) (51) 

for any function / which for functions differentiable at large u gives Asf{u) ~ /'(^) for 
z? > 5. 

This results in 

Oa = ^«(|74)'^' ■\/^(2sgn3/2(^)®id-3(x3/2(^) + ^"'sgn3/2(^))®P 

+ p-\x3/2{^) - '^P~'sgn,/,{P)) ®P' + OiPyu')) (52) 

and 

4 = 1^1 (2A3|^>r/' ® id + 3{P-'A,\P\'/' - A,{\P\-'/hgn{P))) ® P 

- |(0-iA3(|Pr^/'sgn(0)) + iAal^r^'/') ® + ^(PV^')) ■ (53) 



where z/ = 6p/{'yip) and P = Ap^jijlp). From these expressions we can read off the first 
coefficients 0^^j,„ and Ofl^. 

Note that, while the expressions x&i sgn^ and A5 can be directly extended to be applied 
to operators, the expansions are not densely defined on T^pert because of inverse powers of 
u. In fundamental expressions where the classical analog has inverse powers of p, such as 
matter Hamiltonians, one can still define densely defined operators [2^ I24j . Here, however, 
we started with densely defined operators and obtained inverse powers only after expanding 
in the perturbative sector. We thus cannot change those terms and make them densely 
defined. For what follows we can mostly ignore this issue, disregarding states at small v 
on which the operators cannot be applied. We will come back to it in the discussion of 
singularities in Sec. El 

Next we must have the action of the several terms to the left of Oa and Oc in . Up 
to order we have 

H = A7r~^i^~^G'^£p'^ 2 {sin^ |ccos^ |c id + | (- sin Iccos^ + sin^ ^ccos ^c) (g) e 

+1 (—3 sin^ |ccos^ |c — cos'' |c — sin^ ^c) ® £^ + ^ (sin ^ccos^ — sin^ |ccos |c) ® e 
+0(£^)}0a + 

{sin^ iccos^ |c ® id + (sin locos'^ |c — sin^ |ccos |c) ® e 

+ j (—6 sin^ |ccos^ |c + cos'* |c + sin'' ^c) ® + | (— sin ^ccos^ |c + sin^ |ccos |c) ® 
+0{e')]d}. (54) 



3 
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From ()25|) . the relevant combinations of sines and cosines on the isotropic basis above 
have the matrix form 



(z/| sin^ ^ccos^ |c 
{v\ sin iccos^ \c 
sin^ |ccos \c 

/-I • 4 1- 

(z/|sm 

/-I 4 1- 
(Z/| COS -^c 



l/i) = (5p,/i+4 - 45p,^+2 + 65p,^ - 45p_^_2 + 5p,/i-4) , 
= 2"^ (5p,/i+4 + 4(5p,^+2 + 6(5p,^ + 4(5p,/i_2 + 5p,/i-4) • 



(55) 



So the total action in is given by 



mm) ^ 4. 

+ '^P,/I+4 



(4 - 18£^) Oa(z>) + 20c(i; 
(-2 + 2i€ + e^) OAii^) + (-1 - 2i€ + 2e^) 6^ 
(-2 - 2ie + e^) 6a{v) + (-1 + 2ie + 2£=^) ^(z?) + 0[e^ 



(56) 



as matrix elements in ?-^iso taking operator values in Tis (for fixed z/, Oa{i^) and Oc{v) are 
considered as operators on |z/) ® — '^s)- 

We can act with (jSEI) on a generic state |s >= Sy{e)\u > and use (j^Sl ES))- For this 
we write the Hamiltonian in the form 

< z/,e|^|s >= {^+4Sp+4(£) +SpSp(£) +C'p-4Sp-4(£)} , (57) 



where 

Cu = 



-i(l^r/'sgn3/2(z>) + Iz^lAslz^r/^) + l{H''sgn,„{u) - 2\u\A 
+ {lH'\xV2{y) + z>-'sgn3/2(z>)) - sgnz7A3|z>|^/2 ^ Agdz/l-i/^sg^^^)) p 
+ 1 (-|z?r^/^sgnPx3/2(z>) + i|z?r^/^sgn3/2(z?) + sgnz/A3(|z/|-i/2ggj^z/) + P 
- (|z>r/\3/2(z>) + |z>r'/'sgnz7sgn3/2(z>) + 2sgnz/A3|z>|'/' - 2|z/| Agdz^r'/'sgnz/)) leP 
+ 1 {\u\^l\gn,j,{u) + 4|z7| A3|Pr/^) + 0(3) , 



(5^ 



-1/2 



-2(|z?r/'X3/2(z^) + |z^| 

+3 {\v\'^/'^sg-n.vxz/2{^, , 
-6|P|i/2gg^3/2(z>)£' + 0(3) 



sgni/sgn3/2(z^) - sgnz/Aalz/]^/^ + |z/|A3(|z/| ^/^sgnz/)) P 
^/2sgn3/2(z>) - sgnz/AgdPl^^/^sgnz/) - 

(59) 



and Ap = Cp. Note that this is not the adjoint of Cp but the complex conjugate (acting only 
on the numerical coefficients, not on possible factors of the imaginary unit in a derivative 
representation of P). Complex conjugation and adjoint are not the same in this case 
because of the eP-term. Also the whole constraint operator is not self-adjoint since we 
started with a non-symmetric ordering with holonomy operators to the left. 



15 



4 Evolution 



The constraint equation given by the operator (|K7j) . i.e. requiring physical states to be 
annihilated by the constraint operator, results in a difference-differential equation for s„(e) 
which can be interpreted as evolution equation in the volume v as internal time. This 
equation is difference-differential because we keep the underlying discreteness from quan- 
tum geometry in the isotropic variable, but treat the anisotropy perturbatively. Also in 
isotropic models one can have difference-differential equations if there is a matter degree 
of freedom such as a scalar, but the crucial difference here is that there are differential op- 
erators at all levels of the difference equation and in particular the leading ones. This has 
implications for the recurrence scheme since those operators have to be inverted in some 
manner. We perform here only a basic analysis to check if the scheme has satisfactory 
properties concerning evolution and relation to classical behavior. 

We focus on the stability issue for which in our context three different notions occur. 
First, in the perturbative setting stability is already relevant at the classical level where 
one needs to make sure that anisotropies stay small enough for sufficiently long time. As 
can be seen from the discussion leading to ()22j). this can be ensured for some range of 
evolution at least in the vacuum case for our situation. The second notion is stability at 
the quantum level in the sense of requiring that there are no exponentially growing 
solutions to the difference equation. This will be the focus in the following sections. Finally, 
also numerical stability issues can arise when difference or difference-differential equations 
are solved numerically. 

We consider the evolution equation (j37j) without matter for large for which the 
constraint equation reduces to 



23/2 ^3/2 



4 4 
3^3 



_2 6 
3~ V 



2 

3^ 
2 



^p2 

2z/2 
3 



3 3 3 V Iv^ 

1--P- I2e' - , ^^^^^^ 



(60) 



Now we propose a restricted function space for the solutions of the form 



(61) 



where 'p]f' are complex valued polynomials of order k. Since we expanded the Hamiltonian 
constraint to second order, we also should disregard terms of higher order in the solutions 
and work with k = 2. Then, 



(62) 



where Aq, Ai and A2 are complex valued functions of the discrete time u only, to be 
determined from initial values through the difference equation. This class of functions is 
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general enough for our purposes because, even though Eq and a are fixed, the Gaussian can 
be deformed and spread if Ai and A2 grow as functions of u. Moreover, to the given order 
of expansion also non-zero expectation values for are allowed through complex valued 
Ai and A2. For a Gaussian state peaked at anisotropy p° we need a phase factor e^*^^"'''^^? 
which for small enough p° can be expanded, only contributing imaginary parts to Ai and 
A2. For this, we need to assume p° <^ ip which is possible but not guaranteed by our 
general perturbation assumptions. Still, for the tests we are going to perform such small 
values are sufficient. For more detailed aspects one has to generalize the allowed class of 
functions. 

Computing (IHUjl using P = —^id/de and (jUTj) . (jU^ for large z/, and collecting terms of 
the same order in e we obtain 

Ao{u + 4) {2a^ - 8z>V^ - el) - 2Ao{u) {2a^ - Su'^a^ - 2iua^eo - el) (63) 

+Ao{u - 4) (2(7^ - 8z?V^ - el) 
+Ai{u + 4) (-4(7%) - 2Ai{u) {-Aiua^ - 4(t%) + Ai{v - 4) {-Aa'^eQ) 

+A2{p + 4) (-8a^) - 2A2(z7) (-8a^) + A^i^p - 4) (-8a^) = 

from leading order in e, 

Aq{u + 4) {m'^a^ - 12P(T% + 2£o) - 2Ao(z>) {-2iua'^ + 2^0) (64) 

+Aq{v - 4) [-m'^a^ + 12z/a% + 2^0) 
+Ai(z7 + 4) (6(t2 - 2Aua^ - 8i?V^ - - 2Ai{u) (6(7^ - 8z>V'^ - eg _ 2iz/or%) 

+Ai(z/ - 4) (6(t2 + 2Aua^ - 8z/V^ - e^) 
+A2{P + 4) (-8a2£o) - 2A2{iy) (-8zP(t^ - 8(t2£o) + A2{iy - 4) (-8(t2£o) = 

from linear order, and 

Aq{u + 4) (-1 + I2ua^ + Au'^a^) - 2Ao(z/) (-1 + 36z>V^) (65) 

+Ao{u - 4) (-1 - Uua^ + 4z/V^) 
+v4i(z/ + 4) (8iz/V^ + 2eo - 12 - 2Ai(z/) (2iz/cT2 + 2eo) 

+Ai{u - A) {-Siu^a'^ + 2eo + Uua^eo) 
+A2iu + 4) (lOd^ - A8ua^ - 8z>V^ - eg) - 2^2 (z>) (lOa^ - 8z>V^ - 2iua^eo - eg) 

+A2(z>-4) (lO(T^ + 48z>(T'*-8z/V^-eg) = 

from quadratic order. As the equations are quite complicated, we discuss here mainly their 
approximation for large volume which can already be used to see relations to the classical 
behavior. 

4.1 Asymptotic T> solution 

While the perturbation scheme does not require a relation between the magnitudes of u 
and £0 or cr ~ £^0 ^ 1, for a given choice of eo and a consistent with the perturbation 
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assumptions we can look at the asymptotic form of the equations where only highest 
powers of v are used. For the second order equations, this requires v ^ cr^^ which is our 
assumption for the asymptotic regime. Simplifying the equations then leads to 



Ao(z> + 4)-2Ao(z>) + Ao(z>-4) 
Ai(z/ + 4) -2Ai(z/) + Ai(z/-4) 

^2(1^ + 4) -2A2(z>) + v42(z/-4) 



0, 

^[Ao(z> + 4)-Ao(z7-4)], 

i \A^{v + 4) - \%A^{v) + A^{v - 4)] + 

^[Ai(z7 + 4)-Ai(P-4)]. 



(66) 



We can already observe that, in contrast to (jHSl), these equations are independent of the 
value of £0 in complete analogy to our previous observation for the asymptotic classical 
solution (j22|) . Also, the equations have imaginary terms which is necessary for to change 
in internal time (in the ^-representation chosen here, enters through the phase of the 
wave function). Even if p° is initially zero, i.e. Aq, A\ and A^ are all real, the wave 
function must become complex after some steps such that anisotropy also in the triad 
arises automatically. 

Before solving (jUBj) we generalize it to allow dust like matter (matter Hamiltonian 
independent of v) such that the middle coefficients on the left hand sides of the equations 
are 2p with a constant p < 1, reproducing the vaccum case above for p = 1. We then 
put the system in a first order form, considering solutions only on the lattice 4Z which is 
sufficient to see the behavior of long-time evolution. Accordingly, we introduce the vector 



v(z/ + 4) 



With this definition we can cast 
four) : 

M v(z/ + 4) + v(z/) = 
where the matrix Q is 



/ Ao(z^ + 4) \ 
Ai(i/ + 4) 
v42(z/ + 4) 

A,{y) 

\ M^) J 

in a first-order linear matrix equation (of stepsize 
or v(z/ + 4) + g v(z/) = 0, (67) 



Q 



/ -2p 

-2ip -2p 
9 + p 2ip -2p 
-10 
0-10 
\ 0-1 



1 \ 

2i 1 

-2 2i 1 





/ 



It is possible to bring this matrix to its Jordan canonical form (see App. \^ J = S ^QS 
where 
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/ Ai 1 \ 

Ai 1 

Ai 

Aa 1 

As 1 

\0 OA2/ 

Here, Ai,2 = ± y/p^ — 1 — p are the two eigenvalues of Q with muhiphcity three (which 
degenerate further in the vaccum case p = 1). 

From here we can obtain the solution starting from initial conditions in ()67|): 

v(z7o + 4A;) = gV(z7o). (68) 

We obtain 

= SJ''S~\ (69) 
where J is the Jordan form for Q above and each one of its block powers is 

is the off-diagonal term in J. 

We obtain the transition matrix S by computing the generalized eigenvectors associated 
to Ai_2, which we do not write down; however we write the general /c-power of J above: 



( A? 




©At- 











\ 





x\ 





















\\ 





















A^ 




















XI 


kXt' 




V 
















I 



From this the solution for the A functions is 
Ao(z>o + 4fc) = l(p2 - l)-i/2 [{-X\X2 + AiA^)Ao(z>o) + (A^ - )A)A^{u, - 4)) , (71) 
Ai(z7o + Ak) = ^(p2 - l)-i/2 [^^_x\X2 + AiA^) A(z7o) - ik{-XlX^^-^ + X\-^Xl)A^{u,) 

+ i{Xl -X\- kX\~^X2 + A;A^"^Ai)Ao(z7o - 4)(A^ - X\)A^{uo - 4)) (72) 
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and 

^(z7o + 4fc) 



+ 



+ 



+ 



+ 



+ 



+ 



9 



_(p2 _ 1)3/2 



-1 + 18p + 2p2 



v9 



1 



fc\ 2A^ 



2yp2 



k-2 

^ (l - 5p2 + 4p^ - 3pv/p2 - 1 + 4pVp' - 1 



-9 - 2p + 18p2 + 2p3 ^ -1 + 18p + 3p2 + 18p3 - 2p^ 



\ fc-2 
^2 



p2 - 1 

2 - 12p2 + 18p^ 



(p2 - 1)^ 



+ 



ik ( xt^xl- X'^-'Xi 



, fc-1 • 



1 



^i(^^o) 



(p2 - 1)3/2 

6p - 14p3 + 8p5 

(p2 - 1)3/2 
-X1X2 + A2A1 



At 



l_9p + p2 i + 9^_^2 

+ ^2" 



(p2 - 1)3/2 



(p2 - 1)3/2 



+ 3A; -Af 



3p 



+ 



3 + p 



A^ 



fc-2 



+ A 



fc-2 



4p 



(p2-l)f 



4p + 2 
2 

A; 



2p2-l 



+ A 



fc-i 



-1 - 



3p 



+ 



2p^- 



(p+l)2 (p+ 1)3/2^ 

2 — Aj , ^1 ^2 



Ao(z>o - 4) 
Ai(z/o -4) 



2v/;^^ 



A2(z/o -4). 



(73) 



This solution is valid only for |p| 7^ 1, for otherwise the eigenvalues would be completely 
degenerate. 

So we have obtained the solution to the evolution equation in the asymptotic u regime. 
We should note that the solutions grow as powers of the order k (no larger than second) 
as a consequence of degeneracy, although the leading coefficient Aq is bounded. 



4.2 Inclusion of matter 

One can hope that matter terms remove the degeneracy present for the case for pure 
gravity, which similarly happens in the isotropic model (see jTJ for an explicit isotropic 
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vacuum solution which is also unbounded). As the previous discussion showsm however, 
this is not the case for dust which has a constant diagonal Hamiltonian, H = EI. 
As the next simple possibility we can consider a cosmological constant term 

^matter = AV^ . (74) 

Using (jlHj) the evolution equation for large u fl60|) modifies to 

4 4. 2 2 6. ^ 3 

3 3 3 u Zu^ J 

4 4. 2 2 6. ^ 3 
3 3 3 2//"^ / 

_ 2p _ 12^2 - lp2 

- ^/') ■ 

With the same ansatz (jHH), (|S21) the right hand side of each equation becomes 



6 

for (jnni), 

Z73/2 



[Ao(P)(54(T^ + 32z/V^ - 275^) + Ai(z/)(-108(t%) + A2(z/)(-216(t^)] (76) 



6 

for ((nil) and 



[Ao(z>)(54eo) + Ai{u){lQ2a'' + 32z/V^ - 27e^) + A2(z>)(-216a%)] 



6 



[Ao(z>)(-27) + Ai(z/)(54eo) + A2(z>)(270a' + 32z/V^ - 27^^ 



for dnni). 

It is easier to check the degeneracy in the asymptotic v regime: in this case the right 
hand sides of acquire an extra term '^tiG'~^'^ Ku^f^ Ai{v) , respectively. Coefficients in 
the asymptotic difference equation are thus no longer constant and we do not find explicit 
solutions. Still, around any fixed v eigenvalues remain degenerate. 



4.3 Large u behavior 

Rather than looking at more involved matter choices, we now use only large v and keep 
next order terms in v^^ . This gives rise to more complicated correction terms in (j63|l 
compared to the asymptotic equations, changing all coefficients. Rather than trying to 
solve this second order system of difference equations, we perform a local stability analysis 
around large fixed v similar to that in the previous section 



21 



To this end we put the system ()60|) again in a first order form ()67|). where now the 
matrices M, N, Q are functions of {u,eQ,a). In this case the stabihty of the system is 
obtained if and only if the module of the eigenvalues of the matrix Q above are less or 
equal to one. However its characteristic polynomial is of order six and, recalling the form 
of (jnni), each one of its coefficients is a very complex function of {u,eo,(T). Therefore we 
perform a numerical analysis to see the behavior of each one of the coefficient functions 
Ai{iy) in 

The characteristic polynomial is of the form 



which one can solve numerically to find choices of Eq and a where all six roots have an 
absolute value less than or equal to one. 

5 Perturbations and the singularity 

When perturbing in basic variables, the expansion is valid only for -C p and thus 
breaks down very close to the classical singularity. Indeed, the main aim for studying 
perturbations is to understand effects in structure formation where one needs at least 
perturbative inhomogeneities for which anisotropics serve as a model. This is usually 
modelled in semiclassical regimes, in an inflationary phase or sometimes through a bounce 
at relatively large volume, far away from the Planck scale. It is nevertheless instructive 
to discuss implications for the singularity issue itself close to vanishing p, as they might 
appear in a perturbative setting. 

In the asymptotic u equation we saw that the relevant matrix M in fl(j7|l at highest 
order of the difference equation is always invertible such that evolution cannot break down 
in this regime. But this is not guaranteed at small u where a breakdown of evolution 
can occur even at z/ > where one would not expect singularities classically. There are 
even divergences in coefficients of the difference equation such as in Cp for z/ = 3 coming, 
e.g., from A3|z/|~'^/^. This is certainly a small value for u, but note that can be small 
compared to p even for z/ of the order one which perturbatively is consistent using the 
fact that the perturbative operator has a continuous spectrum. As discussed earlier, 
small do not automatically bring us into discrete regimes, unlike small p. One can 
thus set up the perturbative evolution such that all assumptions for the perturbation 
scheme are satisfied, and yet the evolution breaks down at non-zero volume. (Starting 
at large volume, the perturbative approximation is expected to break down much before 
the classical singularity is reached. But generically, a breakdown of evolution cannot be 
ruled out since there are initial conditions for which the evolution breaks down before the 
perturbative scheme ceases to make sense.) 

There is thus no removal of singularities in the perturbative model. Even worse, the 
perturbative equation could make us believe that the non-singular evolution of the isotropic 
model is very special and could not extend to less symmetric models, owing to the fact that 



A*^ + + /i4 + h \^ + h2 X^ + hi \ + ho 








(77) 
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for e = P = all divergent coefficients disappear. In such a situation, it is necessary to 
use the full model from which the perturbative equation is derived in order to see whether 
or not it breaks down. As it turns out, for loop quantum cosmology the evolution is 
non-singular even for anisotropic models, in contrast to the perturbative appearance. 

A different source for a breakdown of evolution appears because coefficients of the differ- 
ence equation are now differential operators on the perturbative degree of freedom. These 
operators must be inverted to proceed with the evolution, which is not generically possible 
and may require special boundary conditions. Note that, if a physical inner product is not 
known, even generic invertibility would not be enough: one needs to make sure that any 
initial state can be evolved without breakdown to conclude singularity freedom. Even if 
almost all initial values would lead to non-singular evolution, all of them could be ruled 
out by the physical inner product (see also P). We do not discuss this in detail here but 
only conclude that, while the fully quantized anisotropic model is non-singular, there is no 
such statement in the perturbative quantization. 

It is interesting to compare this situation with recent results in string theory [26^ where 
inhomogeneities on a background are seen to prevent the occurence of a bounce instead 
of a singularity (which would otherwise be possible in the corresponding homogeneous 
model). This looks similar to our perturbative quantization of anisotropies on an isotropic 
background, even though the models and techniques are certainly very different. In con- 
trast, non-perturbative background independent models studied so far are non-singular, 
including inhomogeneous ones which classically have local physical degrees of freedom [2Zj . 

There is a further disadvantage of using perturbative or semiclassical degrees of free- 
dom in order to discuss the singularity issue, as for instance suggested in |2H]- Such a 
perturbative treatment is unlikely to remain valid close to a singularity where potentially 
all degrees of freedom can be excited strongly. Perturbations certainly allow one to in- 
clude all degrees of freedom which is important for phenomenology, but properties of the 
singularity can be extremely blurred as we have seen here. In contrast, symmetric models 
completely remove many degrees of freedom which may have to be reinserted later on when 
discussing their effect on evolution. But symmetric models themselves are often classically 
singular, and those singularities can be studied by quantizations of models. With the mod- 
els now available, the most characteristic types of classical singularities can be studied also 
at the quantum level. This can then show how classical singularities can be resolved by 
quantization, and there is indeed a general mechanism in loop quantum gravity. 

While the fundamental singularity issue concerning the extension of solutions to the 
Hamiltonian constraint equation cannot be analyzed perturbatively, this is still possible 
for phenomenological models of bounces as studied in j2H] and used for the construction 
of oscillatory universe models in [20]. Such bounces need to happen sufficiently far above 
the Planck scale in order to ensure the validity of effective equations , and thus avoid 
a potential breakdown of the perturbative discrete evolution. This will also be true for 
perturbative analyses of inhomogeneities. Also at the inhomogeneous level, one can employ 
approximate but non-perturbative methods in the spirit of the BKL picture [22] ■ In con- 
trast to perturbative schemes, this can then also be applied to the fundamental singularity 
issue, as done in the loop quantum gravity context with preliminary investigations in . 
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6 Conclusions 



We have mainly illustrated the reduction procedure to symmetric situations available for 
loop quantum gravity, and contrasted it with a new method for perturbative quantum 
degrees of freedom in the case of anisotropy. We consider this mainly as a model situation 
for the cosmologically more interesting but technically much more complicated case of 
perturbative inhomogeneities. This already allowed us to draw some cautionary conclusions 
for the singularity problem when discussed in a perturbative setting, as well as for the 
relation between symmetric models and perturbative less symmetric ones. 

The resulting evolution equations are of a new type compared to other models in loop 
quantum gravity. Although not studied in detail here, they are treatable at least by 
numerical schemes and provide an interesting setting to understand evolution equations in 
general and numerical techniques to solve them. New schemes for an efficient solution still 
need to be developed, also taking into account the issue of pre-classicality (Mj ignored here. 
This condition requires solutions not to vary strongly on the scale of discreteness. It can 
be expected for semiclassical behavior, but is also necessary to ensure the correct behavior 
of wave packets. Techniques to extract pre-classical solutions, such as those based on 
generating functions continued fractions or numerical ones using suitable function 
bases ^3 > have been introduced but are not yet available for the equations considered here. 

Also the equations themselves can be generalized in several directions. Within the same 
degree of symmetry one can discuss different matter choices and, at the level of the ansatz 
(jUT|) . more general classes for the functional behavior of the wave function depending on 
anisotropy. Moreover, different orderings of the constraint operator, most importantly the 
symmetric one, can be tried. The bigger step consists in including inhomogeneous de- 
grees of freedom perturbatively for which the analysis in this paper serves as preparation. 
Classically, the split into symmetric and non-symmetric variables would be given by a 
mode decomposition with respect to the symmetric background. A priori, there are no 
insurmountable difficulties in repeating the analysis here for this situation, even though 
technically it would be much more involved. Inhomogeneities can already be included 
explicitly by using midi-superspace models for which the loop framework has been con- 
structed in ^Sl EH- These situations already allow the inclusion of inhomogeneities of 
cosmo logical interest. 

Our analysis also shows that the reduced constraint of the more symmetric model still 
plays a role for the less symmetric one and is only amended by perturbative correction 
terms. Those terms are derived by expanding the less symmetric constraint operator, 
and not by quantizing the classically expanded constraint in the fashion of j^H]- Indeed, 
the additional terms are much more contrived than expected classically. Although more 
complicated, these terms are derived from the less-symmetric model and are thus well- 
motivated. After a detailed study of the evolution, such as restrictions from quantum 
stability, one can draw conclusions for the form of the less symmetric constraint operator. 
(Such restrictions indeed arise easily as the closed isotropic model, jTU] compared to [iO] . 
and the more detailed analysis in [23] showed.) Eventually, with a better understanding 
of the relation between models and the full theory also at the level of the Hamiltonian 
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constraint, one can then hope to restrict possible full constraint operators and reduce their 
ambiguities. 
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Appendix 

A Approximate action of operators on peaked states 

For the perturbative expansion of the Hamiltonian constraint at the quantum level it is 
necessary to define and discuss an appropriate notion of expansions of operators. We are not 
interested in expansions valid for any state of the Hilbert space but only in expansions that 
can be used for states corresponding to small anisotropy such as those used in (jHH). Thus, 
the expansion parameter is not contained directly in the original operator but introduced 
through a restricted class of states considered for the action. We then start by choosing a 
space Vs of states in the original Hilbert space H, or even a subset of a cylindrical space 
Cylg, depending on the intended perturbation parameter S (or several ones such as Eq and 
p° in the main part of this paper, and define 

Definition 1 An operator Ok on Ti is an approximation of order k to an operator O if 
Ok^ agrees, with respect to the inner product, with for all ^ E T>s up to order 5^^^^ , 
i.e. 

(vi/'|(6-6fc)vi/) = o(5'=+i) 

for all G Cylg. Agreement with respect to the inner product in this sense will he denoted 
by Ok^ ~ 6\&. 

To show that this is realized for the expansions performed for anisotropics in this 
paper, we choose Vq^^pO to be the span of states of the form = cr~^/^s(g) exp(— (g — 
go)^/4o"^) exp(ip°g) for a given a < qo and a polynomial s (we do not normalize states 
but keep track of factors of a). Expansions performed here always lead to approximations 
Ok which are polynomials in the basic operators q and p of the anisotropy. This implies 
that any Ok considered here maps the space to itself and can be used to show that the 
validity of expansions is preserved by taking products: 
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Lemma 1 Let O^^ and O^^ be approximations of O^^^ and 0^'^\ respectively, such that 
they fix the spaces T>s and Cylg. Then O^^O^^ is an approximation of O^^^O^^K 

Proof: By assumption we have 6(1)6(2)^ ^ d(i)6f b ecause we can simply replace 
by in the definition of agreement. Since is in Vg for G we have 

which finishes the proof. □ 

We can therefore split the expansion into different steps as done in the main calculations 
here. This also simplifies the proof that the performed expansion of the Hamiltonian 
constraint is an approximation because we only need to show that multiplication operators 
can be approximated by simply expanding the multiplying function, which we would apply 
for our calculations in both the coordinate and momentum representation. We then only 
need 

Lemma 2 Let O^'(g) = 0{q)'^{q) be a multiplication operator on H = L^(]R, dg). Any 
partial Taylor sum Ok{q) — Ylt=o S'^^s* '^^^c ^0 an approximation Ok'ifiq) — 

Ok{q)'^{q) of O order k with respect to Vs as above. 

Proof: Let ^'(g) = cr"^/^s(g) exp(-(g - qof/Aa^) and G Cylg. We have 

/■go+i 
J go 

POO 

+2C2(7-V2 / exp(-(5 - q^f/Aa'')q^+^dq 

Jan+l 



dq^'{q){0{q) - Ok{q))s{q) exp(-(g - qo)VAa^] 



with positive constants Ci and C2 depending only on the supremum of the real and imag- 
inary parts of ^'{q)s{q). The first integral can be estimated by 



rgi 
J go 



90+1 pi 

exp(-(g-go)V4(T')?'^'dg = / exp(-gV4a2)(? + ?o)'+Mg 



a^^^ / exp(-gV4)(g + go/(T)'+Mg 

poo 

< exp(-gV4)(g + go/a)'^+Mg 

^0 



where the last integral exists and does not depend strongly on qQ and a since they are 
assumed to be of the same order. The whole expression is thus of order a^^'^. For the 
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second integral, we have 

/■oo pco 

/ exp(-(g-go)V4fT')?'"''dg = e'^^-rW exp(-(g2 + 2g)/4a2)(g + go + 

Jqo + l Jo 

POO 

< e-(2-)-' / exp(-gV4fT')(g + 2)'=+Mg 

where the last integral is a polynomial in a and the whole expression is thus of the order 
g-(2o-) times a polynomial in a. By our assumptions on the states in T)^, we have a of 
the order or smaller than go which concludes the proof. □ 



B Jordan canonical form 



B.l Generalized eigenvectors and transition matrix 

In this section we briefly review the Jordan canonical form of a square complex matrix. 
For further details and the proofs of the theorems see for example jH]. We begin stating 
the following 

Theorem 1 If a square matrix A of order n has s linearly independent eigenvectors, then 
it is similar to a matrix J of the following form, called a Jordan canonical form 



/ Ji \ 

J2 



J 



J = S-^AS 

, Js 

in which each Ji, called the Jordan block, is a triangular matrix of the form 

( K I \ 



(75 



Ji 



\ 



1 

A. / 



(79) 



where Aj is a single eigenvalue of A. Each Jordan block corresponds to a linearly indepen- 
dent eigenvector. 

This theorem states that we can bring any square matrix to its Jordan form. The 
usefulness of this form lies in the fact that it is easier to obtain powers of the matrix A, as 
can be seen by taking a look at the form of J in (f75j) . (f?^ . The matrix 5* in (fTHj) is called 
the transition matrix. Now we state the procedure to obtain both matrices J and S. 

The basic elements to obtain the Jordan form are the generalized eigenvectors, replacing 
ordinary eigenvectors of a diagonalizable matrix. 
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Definition 2 A nonzero vector x is said to be a generalized eigenvector of A of rank k 
belonging to an eigenvalue A if 

(A-A/)^x = and {A-Xlf'^^y^O. (80) 

If A; = 1 this is the usual definition of an eigenvector. For a generalized eigenvector x of 
rank larger than one belonging to the eigenvalue A define 

Xfc = X, 

Xfc_i = (A-A/)x = (A-A/)xfc, 
Xfc-2 = (A-A/)2x = (A-A/)xfc_i, 



X2 = (A-Xlf-^x = (A-A/)X3, 
xi = {A-Xlf-^^ = (A-A/)x2. 



^11 



One can check that each one of the x^ above is a generalized eigenvector of the eigenvalue 
A. 

Definition 3 The set of vectors {xi, . . . ,Xfc} is called a chain of generalized eigenvectors 
belonging to the eigenvalue A. 

Now the transition matrix S can be constructed from the chains of linearly independent 
generalized eigenvectors of A 



B.2 Powers of square matrices 

The Jordan canonical form of any square matrix A enables us to compute its powers A^. 
Let J be the Jordan canonical form of an arbitrary n x n matrix A such that 

/ Ji \ 

S-^AS = J = 



Js J 



where S is made of the generalized eigenvectors of A and Jg are Jordan blocks. 
Since we have 



\ 



for k = 1,2,..., it is enough to compute for each Jordan block J. An m x m Jordan 
block J belonging to an eigenvalue A of A may be written as 

/Al 0\ /01---0\ 

■•• ■•• 

: 1 



J 







A 1 

A y 



AL 



+ 



\0 



Al + iV, 



/ 
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where I is the identity matrix and is nilpotent. Clearly we have 



(82) 



j=0 



But iV^ = for /c > m and thus, by defining = if /c < / 



J' 



k 




( 



m — 1 



k 



) 



(83) 
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